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ABSTRACT : 

A flow configuration of an existing low-speed, axial flow compressor 
test facility at the Turbopropulsion Laboratory, Naval Postgraduate 
School, Monterey, California, was redesigned in order to improve com- 
pressor inlet velocity distribution, improve mass flow measurement, and 
increase the power input. The design of the new configuration, consist- 
ing of an inlet nozzle, ducting, throttling device upstream of the com- 
pressor, and a diffusor downstream of it, is reported. The redesign is 
based on using the existing free-vortex type and solid -body type bladings, 
and the newly designed symmetrical blading. 
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INTRODUCTION 



The three-stage, axial flow compressor test rig of the Turbomachinery 
Laboratory of the Department of Aeronautics, Naval Postgraduate School, 
Monterey, California, was acquired from the California Institute of 
Technology, Pasadena, California* It was used in the acquired config- 
uration for several test programs until a recent experimental program 
to determine the influence of blade loading on tip clearance effects 
was begun which required the redesign of the facility. 

The object of this report is to present some of the reasons that led 
to the decision to redesign, and to present a review of the design 
criter ia. 

DECISION TO REDESIGN 

The original configuration presented several limitations. See Figure 
1. In particular, the flow rate was controlled by a throttle at the 
compressor discharge, which caused asymmetric velocity distributions 
leaving the compressor. Asymmetric velocity profiles were also de- 
tected at the entrance portion of the inlet duct. Since mass flow 
rate is determined by the integration of the entrance velocity profile, 
this results in an uncertainty in the mass flow rate. Also, inter- 
stage measurements are restricted to a small angle of the periphery. 

A uniform velocity profile ahead of the compressor is required to 
assume axial symmetry and to apply these measurements to the entire 
periphery for overall performance determinations. In the investiga- 
tion of tip clearance losses, an accurate determination of the mass 
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flow rate is essential, and therefore an alternate to the original 
design was required. 

DESIGN CRITERIA 

The facility design was changed to a system in which the mass flow 
rate is controlled by a throttling device located in the inlet duct 
to the compressor. The overall requirements of the inlet section 
are to 

(a) provide an axisyrametric flow as uniform as possible 
ahead of the compressor 

(b) provide the means for an accurate determination of the 
mass flow rate 

(c) allow an easy adjustment of the mass flow rate 

(d) minimize pressure losses at the zero throttling posi- 
tion so that as high an operating mass flow rate range 
as possible may be achieved. 

Requirements (b) and (d) are contradictory. Using an inlet nozzle 
for the flow measurement requires sufficiently large pressure drops 
for accuracy, whereas as low a pressure drop as possible is desired 
to obtain the maximum flow rates. 

The component arrangement is shown in Figure 2. The required inlet 
sections include the inlet nozzle, the throttling device and the re- 
quired ducting between the inlet nozzle and throttle and between the 
throttle and the compressor. The change in design also includes the 
addition of a diffusor section immediately downstream of the compres- 
sor discharging to atmosphere, a new drive motor with increased 
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horsepower capability, and a torque meter on the compressor drive 
shaft. 

With higher horsepower available, higher pressure ratios are attain- 
able and hence an improvement in the measuring accuracy is possible. 
In the following, the design and specification of the components 
cited above will be discussed in detail. 



OPERATING RANGE 



The operating range of the compressor was determined by Vavra (1) 
in the design of the symmetrical blading to be used in one, two, and 
three stages. The compressor speed for each of these configurations 
is based on the absorption of 135 hp at the design point. The drive 
motor is rated at 150 hp, 1180 rpm, 44Ov/30, 190 amps continuous 
operation, which allows for operation at overload conditions. A 
fixed V-belt pulley drive is used to obtain the desired compressor 
speeds of 2290, 1818, and 1588 rpm (1). 

The off-design performance of the sym metrical blading was then de- 
termined. Several simplifying assumptions were made: 

1. The performance of the stage element at the blade mean 
radius, R^ is representative of the whole stage. 

2. The axial components of the absolute velocity ahead of 
and after the rotor blades at the mean radius, R^, are 
set equal to each other with a magnitude equal 

to the average of the original components, i.e., 

V a * 1/2 ( V al + V a2 >R ' 

m 
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3. The energy input to the rotor stage at the design point, 
as indicated by the magnitude of the change in the tan- 
gential component of the absolute velocity, AV , is set 
equal to the value obtained by Vavra (1). 

4. The radial shift of the streamlines can be ignored. 

5. The stator discharge angle remains constant at a given 
compressor rotational speed for all mass flow rates 
within the range of interest. 

6. Changes in the rotor discharge angle are calculated 
assuming the two-dimensional cascade data of reference 
(2) Chapter VI. 

7. The losses in the rotor and stator are equal. 

8. The flow conditions established for a single stage are 
repeated for multi-stage operation. 

The symmetry of the blading, together with assumptions 2 and 3 above, 
modifies the original velocity triangle from that of Figure 3(a) to 
that of Figure 3(b). Note that magnitudes of all angles and veloc- 
ity vectors except Av^ and will be changed slightly during this 

operation. 

A summary of corresponding values at the design point from the orig- 
inal computations of Vavra (1) and the modified values resulting 
from the above assumptions is given in Table I. 

The procedure for obtaining the modified velocity triangle from the 
original design point velocity triangle is given in Appendix A. 



WORK REQUIRED AT OFF-DESIGN CONDITIONS 



The off-design conditions investigated are those operating points at 
a specified compressor rotational speed with differing mass flow rates. 
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Symbolically, these conditions will be represented by affixing the 
superscript prime to the variable under consideration, e.g., v' . 

SL 

Physically, the mass flow rate will be varied by throttling. Analyt- 
ically, the mass flow rate is varied by changes in the axial compo- 

» 

nent of the absolute velocity, V . From Figure 4, note that with 

a 

a fixed direction of the absolute inlet velocity, , that is, a 
constant inlet guide vane discharge angle, , (see assumption 5 
above), and a fixed rotational speed, a) , changes in the axial com- 

i 

ponent, V , are reflected in changes in the direction of the rela- 
a 

tive inlet velocity, , that is, in the angle B^ . It is 
somewhat less convenient to vary B-^ in the analysis than the 
axial component, , so that this means will be used. The inter- 
relationship between the significant geometric and flow angles is 
shown in Figure 5. In particular, note that for a given blade 
geometry, that is, with a given camber angle, 0 , and a given 
stagger angle, y , the relative flow angles B^ and $2 are directly 
related to the incidence angle, i’ , and the deviation angle, 6 ? , 
respectively. Thus, by selecting a b| for an off-design point 
analysis, the incidence angle, i , is readily computed. 

The sensitivity of the deviation angle to changes in the incidence 
angle has been reported by Lieblein (2) for experimental flow in two- 
dimensional cascades. In the flow regime of minimum loss, this 
dependence is given by 



6 1 



6 ref + ( 1 




EQN (1) 
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where the slope, / — J is a function of the solidity, a , and 

the relative flow angle, 0^ . The reference values chosen are the 
design point values from Vavra (1), which were selected for mini- 
mum loss based on this same reference. 

With this value of the deviation angle, the flow deflection angle 
becomes (see Figure 5) 

Ag' = - 6^ = <f> + i' - EQN (2) 



This value of the flow deflection angle is related to the original 

value calculated at the design point. The flow deflection angle 

for the modified velocity triangle at this off-design condition, 

A0 ? , , is assumed to be in the same ratio to the flow deflection 

mod 

angle obtained from the original values, A0 f , as the ratio of the 
flow deflection angle for the modified velocity triangle at the 
design point, A0^ o ^ , is to the flow deflection angle, A0 , from 
the original design point data (1), i.e., 



A0 



ao 1 ao 1 ^modified 

modified A0 



EQN (3) 



Thus, the relative flow direction at the rotor exit 



A6 2 = 0 f 

modified P 1 



- A0 



modified 



EQN (4) 



modified 

The non-dimensional form of the velocities is obtained by refer- 
ring them to the rotor peripheral speed at the mean radius, 



U = R 



m 



These will be designated by the superscript asterisk, 



e.g. , 




a 



V 

a 

U 



V 

a 



R • 

m 
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The modified velocity triangle is shown in Figure 6 . 

At large positive incidence angles, the flow separation on the 
suction side of the blades causes the deviation angles, 6 , and 
the relative flow deflection angles, A 3 , to be significantly 
different than those values calculated above in Equations (3) 
and (4). Howell (3) indicates the magnitude of this difference 
and it has been incorporated into the off-design calculations 
at incidence angles greater than + 2 ° (corresponding to relative 
inlet flow angles, 3^ > greater than 45° •) 

A plot of the resulting calculated relative turning angle ratio, 

777 as a function of the dimensionless incidence angle, 

Ap Ap 

is shown in Figure 7. 

From Figure 6 , note that 



V 



* » 



a 



1 

tan + tan 3 | 



EQN (5) 



and that the increase in the tangential component of the 

JL f 

absolute velocity, A , is given by 

A V*' = V*o - V*i = 1 - V* (tan « + tan g') EQN (6) 

u u z u all 



Substituting the expression for V*' from EQN (5) into EQN ( 6 ) 



results in 



A V* f 
u 



tan 3 ^ 
tan 



- tan 6 ^ 
+ tan g| 



EQN (7) 



The losses attributable to operation at off-design incidence 
angles and the range of operation as defined by the stall limits 
has been estimated using the results of Howell (3) and Lieblein (2). 
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Lieblein found that an equivalent diffusion factor, D 1 can be 

eq 

related to the off-design flow angles as follows: 



D 1 

eq 



cos $' 2 
cos 



[ 



1.12 + a(i' - i) 1 * 43 + 0.61 



cos 



i 

— (tan 3| - tan Bp 



where a = 0.007 for the C-4 circular arc blades. He found that blade 
stall can occur at positive incidence angles for which the equivalent 
diffusion factor, D^ 1 > is greater than 2.0. Therefore, using this 

criterion, the positive incidence angle where stall is likely to occur, 

i f - i 
s t 

ist 9 Can determined . It was found that the ratio, — — , is 

near 0.45. This is also the region where the turning angle, Ag , 
departs from its lower incidence angle behavior. (See Figure 7.) 

Using this value for the positive stall limit together with the refer- 
ence incidence angle established a half range of operation which then 

i * — i 

sets the negative stall limit at a value st = - 0*45 . 

AB 

Howell (3) found that the magnitude of the loss coefficient near stall 
is twice the minimum value. Combining these two criteria for stall 
establishes the magnitudes of the loss coefficients and incidence 
angles at stall. Further, it assumes that the variation of the loss 
coefficient with incidence angle is similar to the observed values by 
Lieblein (2) . This variation, in normalized form, is shown in Figure 
8, where Y is the minimum loss coefficient, and Y f the off-design 
loss coefficient. 



PRESSURE RISE ACROSS A STAGE 



The total pressure loss in a rotor is customarily defined as a fraction 



of the inlet relative velocity head, that is, 

2 



A P t = YJsp w 



m 



and similarly for a stator 

A P = Y h P V? 
t s S in 

where Y and Y are defined as the loss coefficients for the rotor 
K b 

and stator respectively. It is assumed for this analysis that the 
rotor and stator loss coefficients are equal at the design point, 
that is, Y R = Y g = Y. 

For the separate discussion of the compressor, establish the position 
locations as follows: 

Section 1 : rotor inlet 

Section 2 : rotor exit and stator inlet 

Section 3 : stator exit, 

as indicated in Figure 9. 

The overall stage efficiency, ri^ , defined as the ratio of the total 



pressure rise to the theoretical maximum total pressure rise, becomes 

EQN (9) 



p t - p t 
c 3 1 



n t (P - P_ ) 



'1 theo. 



p t - p t 
3 C 1 

p R w(AV ) 
m u 



*2 



At the design-point p = 1 - 

and at off-design p' = 1 - 



— / 

Y ( W;L + V 2 ) 
2x 

Y'(w 1 + V* 

2 x 7 



EQN (9a) 
EQN (9b) 



Combining the dimensionless pressure rise, tt , defined as 

- - C 3 

P i - p * 
t t 

- , 3 1 

ir - = p 0)2 r2 



m 



EQN (10) 
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with the definition of the stage efficiency yields 



- i 



t 



x 



EQN (11) 



3 



Therefore, from Equations (9b), (10), and (11), the average 
total pressure at the stage exit can be obtained in terms of 
the average stage inlet total pressure. 

For a multistage compressor, it is assumed that conditions re- 
peat from stage to stage, such that for Z stages, the overall 
total pressure rise is given by 



Using the above equations, the total pressure ratio, tt 1 , at 

t 3 

the design point is slightly greater than that obtained by Vavra 



the blade height and obviously a better solution than the mean value 
obtained above by using the modified velocity triangle. Thus, the 
total pressure ratio at the design point obtained by Vavra (1) is 
used. Assuming that the off-design total pressure ratios calculated by 
the above modified velocity triangle method can be improved by stating 
the ratio of improved to calculated total pressure ratios at off-design 
is equal to the ratio of that obtained by Vavra to the calculated value 
at design. Thus, the improved (or corrected) off-design total pressure 



t r 



t 




EQN (12) 



(1) , tt^ , which is an integrated average of the variation across 



r 



ratio becomes 




7 T 



t 



3 . 



mod. tt 

t 3 

mod. 



EQN (13) 



imp. 



design pt. 
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For incompressible flow, the off-design power input, L 1 , is related 

to the power input at the design point, L , at the same rotational 

speed by t f 

L ' V* X 

~ " v*~x E Q N 

a 

The corresponding calculations are indicated in TABLE 2. and in 
Figure 10. 

The expected mass flow rates range can be determined from the design 
point values given by Vavra ( 1 ) and the mass flow functions at the 
assumed stall limits interpolated from Table 2. A summary of the 
results is tabulated below: 







i 


2 


3 






Stage 


Stages 


Stages 


Speed 


RPM 


2290. 


1818. 


1598. 


Design Volume Flow Rate 


3 

ft /sec 


1034. 


820. 


718. 


Design Mass Flow Rate 


slugs/sec 


2.452 


1.944 


1.708 


Maximum Mass Flow Rate 


slugs/sec 


3.059 


2.426 


2.131 


Minimum Mass Flow Rate 


slugs/ sec 


1.943 


1.541 


1.353 



FREE-VORTEX AND SOLID-BODY BLADING 

The performance of the free-vortex type and the solid-body type blading 
has been measured in previous tests. The expected operating range 
using these bladings at the desired compressor speeds for the symmetri- 
cal blading was checked for compatibility with the redesigned facility. 
Performance maps of the free-vortex and solid-body bladings are present- 
ed in Figure 5 of Reference (4) and Figure 46 in Part 2 of Reference (5). 
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The rotor tip speed, coR^ , was used in the correlating dimension- 
less parameters, ip , \p' , and <J> . The data presented are put into 
the form used in this report by changing the reference speed from 
the rotor tip speed to the speed at the mean blade radius, R 

m 

In addition, the symbology has been changed to conform with this 
report. The necessary conversions are summarized as follows: 



Flow coefficients: 



Work coefficients: 



X = 



II 


V 

a 


a 


(jOR 

m 


__ 




AH 


.2 =4 



(wR, Y 

m 



m 



R 

m 



EQN (15) 



EQN (16) 



where AH is the increase in the total enthalpy per stage. 

, „ v 2 



Head coefficients: tt 



P - 
to t- 
— ^ 

p(uR m ) 2 ^ 2 

m 



R 



m 



EQN (17) 



Also, the average total-to-total stage efficiencies are related by 



~ \ l 

n t X ip 



EQN (18) 



Figure 11 presents the off-design performance of three stages of 
free-vortex type blading. Included on the figure is the pressure 
ratio of one stage of solid-body type blading. The behavior of the 
solid-body type blading is very nearly that of the free-vortex type 
blading, thus the off-design performance of the free-vortex type 
blading will be assumed to represent the behavior of both types. 

The horsepower, L , required at the design speeds are first determined 
to compare with the available horsepower of the new drive system. 

L = * (R t " R h ) P f Z *t 3 (wR m ) 3 EQN (19) 
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For the redesigned facility, R t - 18.0 in., R h ~ 10.6 in., and 
R m = 14.4 in. Also, jp = 0.5625 ( 0 = 0.45 ) at design conditions, 
and p = 2.731 x 10 ^ slugs/ft^. From Figure 11, 



for 


z 


= 1 : 


7 = 0.258 

3 


and 


n t = 0.897; 


for 


z 


= 3 : 


7 = 0.258 


and 


¥ t = 0.85. 



3 



Assume for two stage operation that tt = 0.258 and that n is the 

3 t 

arithmetic average of n t for one and three stages, i.e., 

T\ = 1/2 ( 0.897 + 0.85 ) = 0.873 

The power required by the free-vortex type and the solid-body type 
bladings at the three rotational speeds fixed by the symmetrical 
blading design for one, two, and three stages operation is summarized 
in the table below: 

DESIGN POINT HORSEPOWER REQUIRED BY THE FREE-VORTEX 
AND SOLID-BODY TYPE BLADINGS 



Stages 

z 




\ 


N 


Horsepower, 
= 1588 rpm, 1818 rpm. 


L 

, 2290 rpm 


1 


3 

0.258 


0.897 




25.1 


37.6 


75.2 


2 


0.258 


0.873 




51.5 


77.3 


154.5 


3 


0.258 


0.850 




79.4 


119.1 


238.1 


From the 


above table, 


operation at 


2290 


rpm with two or three 


stages 



is not possible with the 150 hp motor. In addition, a conservative 
extrapolation for three stage operation at the minimum flow stall limit 
( <p = .30 ) is obtained from that of one stage. This off-design 
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horsepower is approximately 1.8 times the design horsepower. Thus, 
for three stage operation at 1818 rpm, the off-design required 

t 

horsepower, L f = L x — = 119.1 x 1.8 = 214.4 hp. This precludes 

operation at this off-design limit. 

The expected mass flow rates range can be deduced from the data in 
Reference (1). The flow function of Reference (2), 0 , is 

0 - 7 eqn ( 20 ) 

“r R t * < R t - *h > 

The flow function, , of this report is related to 0 as follows: 

- R t 

<P = 0 jp EQN (21) 

m 

A summary of the results is tabulated below: 

MASS FLOW RATES OF THE FREE-VORTEX AND 
SOLID BODY TYPE BLADINGS 





Stages 




1 




2 




Speed 


rpm 


2290 


1818 


1588 


1818 


1588 


Design Volume 
Flow Rate 


3 

ft /sec 


732.3 


581. A 


56A.2 


581. A 


507.8 


Design Mass 
Flow Rate 


slugs/sec 


1.736 


1.378 


1 . 20 A 


1.378 


1.20A 


Minimum Mass 
Flow Rate 


slugs/sec 


0.965 


0.766 


0.669 


0 . 88 A* 


0.7A8 


Maximum Mass 
Flow Rate 


slugs/sec 


1.929 


1.532 


1.3378 


1.720* 


1.405' 


Results for two 


stages are 


estimated as the 


average 





of the flow rates at one and three stages. 
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INLET NOZZLE 



Two requirements of the inlet nozzle are to provide for a suffi- 
ciently large pressure drop from the atmospheric pressure, P , 

a 

to the nozzle throat static pressure, p| , for reasonable flow 
measurement accuracy, along with a minimum distortion of the 
velocity profile. 

For the isentropic flow of a perfect gas with constant specific 
heats, the mass flow rate through the inlet nozzle is given by 

. 3 * 



"is. - * ^ t p a • 2 (P a - P l> ]h l 1 - 47 



EQN (22) 



where the higher order terms of the nozzle pressure ratio ex- 



P - P 

pansion has been neglected. For a pressure ratio, a 1 



= 0.06 



a 

(corresponding to a pressure drop, P - P* , of approximately 25.0 

a x 

inches of water), results in an error of less than 0.05%. 

The pressure drop for the expected minimum flow rate was calculated 

using Equation (22) above with P = 14.696 psia, T = 59°F, y = 1.40, 

a a 

and setting . The minimum mass flow rate of approximately 

0.67 slugs/second occurs for one stage operation of the free-vortex 
type blading at 1588 rpm. The resulting pressure drop of less than 
0.35 inches of water is too small a pressure drop for accurate 
measurements of mass flow rate. In order to increase this pressure 
drop, the nozzle throat area must be reduced. The variation of 
pressure drop with mass flow rate for various throat diameters was 
calculated and the results shown in Figure 12. Also indicated on 
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the figure are the expected mass flow rate ranges of the symmetri- 
cal, free-vortex and solid-body type bladings. Selection of the 



design radius ratio, 




depends on the criterion set for a 



minimum acceptable pressure drop. If a pressure drop of 2.0 inches 
of water is selected as a minimum, the largest 3 permissible for 
the minimum mass flow rate of approximately 0.67 slugs/second is 
approximately 0.65. At the minimum mass flow rate with two stages 
of free-vortex or solid-body type blading, a 6 of 0.7 is required. 

For three stage operation, the pressure drop is slightly greater 
than the desired 2.0 inches of water at minimum mass flow if a 
3 of 0.7 is chosen. However, to restrict the number of inlet 
nozzles to be made, only one 3 less than unity is desired. Use 
of 3 = 0.7 for one stage of operation is satisfactory for most of 
the flow rate range. The compromise made is that there is a rela- 
tively large pressure drop of approximately 12.5 inches of water 
at the maximum mass flow rate. Whether this is acceptable or not 
remains to be determined in the overall pressure compatibility of 
the system. 

For the symmetric type blading, the minimum mass flow rate obtained 
during three stage operation requires a 3 of 0.9 or less. Hence, 
a 3 = 0.7 will be acceptable if the higher pressure drop is compat- 
ible. For operation with one or two stages, a 3 = 1-0 is satisfactory. 
Therefore, the design criteria for nozzle sizes are as tabulated below: 
Operational Mode Operational Mode 

Free-Vortex or Solid— Body Type Blading 3 Symmetric Type Blading 3 

One Stage 0.7 

Two Stages 0.7 

Three Stages 
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0.7 



One Stage 1.0 

Two Stages 1.0 

Three Stages 0.7 



These selections are to be checked for overall pressure compatibil- 
ity. 

An inlet nozzle contour with a gradually decreasing curvature in 
the flow direction, decreasing to zero curvature at the throat of 
the inlet nozzle, is desired. A curve having this property is 
r = a 0 k y where r is the magnitude of a vector r with fixed origin 
0 whose direction is specified by the angle 0 from a given axis 
(Vavra (6), page 297). See Figure 13. The constant, a , determines 
the scale, and the constant, b , specifies the ratio of QP /QP q for 
a particular 6* and a prescribed curvature at P q . The constant, b , 
cannot be determined explicitly and therefore must be found by 
interation. Selecting an outer diameter of the larger nozzle 
( 0 = 1.0 ) twice that of the throat diameter (in this case, the 
tube or duct diameter) ; and an outer diameter equal to the tube 
diameter for the smaller nozzle ( 3 = 0.7 ) , the following results 
are obtained : 





6 


Throat 

Diameter 

inches 


Large 

Diameter 

(AP*) 

inches 


Axial 

Length 

(QP ) 

inches 


QP* 

Q P o 


b 


Small Nozzle 


0.7 


25.200 


36.000 


12.500 


0.4320 


0.6975 


Large Nozzle 


1.0 


36.000 


60.000 


24.000 


0.5000 


0.5920 



INLET NOZZLE PRESSURE DROP 

The ideal pressure drop from the atmosphere at P to the throat pres- 

a 

sure, p| , is given by 

, • 2 

P a - Pl = 1/2 p 



17 



A convenient reference dynamic head, > can be defined in terms 

of the reference velocity, m R . Recall that the average axial 

m 

velocity, V , when referred to this reference velocity, is the 
a 

V 

flow coefficient, (p = — ~ = V* . Therefore, let 

J a) R a 

m 

2 V 

^ref ^ p re f ^^m^ where p = p . Thus, q = 1/2 p (—jr) ^ . 

ref a H ref a f 



Also, any velocity head, , may be expressed in terms of q re £ : 



q i q ref 



i \ z 

3— 1 (p where p_^ = * 



• 2 2 

From the continuity equation: m = p^A^ V ^ = p & tt ( R fc - R^ ) V & 



If A. = tt R. , then q . = q r 
1 1 9 n ref 



2 2 
(R t - R h ) 

R. 2 

1 



p • Note also by 



R 2 

^h 2 2 2 2 

definition, ct = ~ , and therefore, (R t - R^) = R^ (1 - 01 ) 



R t 2^2 

Hence, q ± = q ref — ( 1 - a ) (p 

R # 

A t R t -2 

Noting that tt = ~~n ~ B , q-. = 

A 1 R^» 



EQN ( 23 ) 



FT (1 " “ 2) q ref f 



and thus, (P a - P x ). deal = ^4 < 1 " “ ) <l re f f 
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The actual pressure drop is somewhat greater than this for a given 
flow rate because of the frictional resistance, and is expressed as 



where the loss coefficient, , varies from 0.03 to 0.05 depending 
upon the area ratio. A mean value of 0.04 was selected. 

SUDDEN EXPANSION DOWNSTREAM OF INLET NOZZLE 

Downstream of the smaller throat diameter inlet nozzle, the flow 
expands suddenly to the duct tube diameter with a consequent loss 
in total pressure. Assuming the pressure is constant across the 
cross-section, the pressure recovery at a cross-section downstream 
in the tube where the flow has re-established a uniform velocity 
profile is obtained by applying the momentum equation: 




factual 




EQN (24) 



- (1 - £ n ) — r (1 - a 2 ) 2 q ref fP 
6 




and the continuity equation: 




This yields a theoretical static pressure recovery, A P £ 



th 



A P 



E 



th 
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R i 2 

Recalling 3 = pT" and q_^ = 1/2 , and assuming incompressible 

flow, 

A P =2 q' g 2 (1 - g 2 ) EQN (25) 

th 

or, in terms of the reference dynamic head 

2 

A P = 2 ir (1 - g 2 ) (1 - a 2 ) q , 0 2 EQN (26) 

th g rer J 

Experiments by Souran (7) indicate that the actual pressure recovery, 
including a tube of a certain length, is somewhat less than the 
theoretical given by 

AP E - ( 1 - C E ) 4P Eth 

where E, varies from 0.03 to 0.05, depending on the area ratio and 

£j 

the length of the tube. The experiments include the frictional pres- 
sure drop with the pressure recovery, and therefore the coefficient, 

, includes both effects. Incorporating the frictional pressure 
drop, A equivalent to P £ - , with the sudden expansion yields 

2 

A P E - AP T1 = P^ - P| = 2 (1 - ? E ) ( -| - 1 ) (1 - a 2 ) q rgf f 2 EQN (27) 

3 

Two important experimental results of Souran (7) are: first, an 

almost uniform velocity profile is attained in the tube downstream 
of the sudden expansion (for an area ratio of 0.5 (3 ^ 0.7), this 
takes place within about six tube diameters); second, the maximum 
pressure recovery for an area ratio of 0.5 occurs after about four 
to five tube diameters. The first result is most important when 
zero throttling occurs, since throttling as proposed should aid in 
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the re-establishment of a uniform velocity profile. A uniform 
velocity profile is required for the input to the compressor and 
is more important than the maximum pressure recovery possible. 

Hence, it is concluded that a total length of six tube diameters 
from the exit of the inlet nozzle to the compressor inlet is 
desired. The division between upstream and downstream of the 
throttle will be determined after examining the throttle charac- 
teristics . 

DOWNSTREAM DUCT LOSS 

The pressure drop in the duct tube between the throttle and the 
compressor inlet is given by 

L 

AP - p 1 - p 1 - r n — 

r T2 3 r 4 ^T2 q T2 2 R t 

where £ ^ is t ^ ie applicable friction factor. Using Equation (23), 

^*2 2 N ^ 2 

A P T2 = £ t2 ( 1 - a ) q ref y E Q N (28) 

INLET TO COMPRESSOR 

The ideal pressure drop associated with the contraction of the flow 

2 

area from that of the tube, tt R , to the annular area of the com- 
2 2 

pressor, tt (R fc - R^) , is given by 

‘V, - (! 4 - P 5> ideal = «5 ' "4 

ideal 
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t 2 | 2 

From Equation (23), = q rgf ; = (1 - a ) q rgf <j) 

Thus , 



A P 



i - a - u-« 2 ) ) , re£ j0 2 

ideal 



The actual pressure drop is greater than the ideal, as represented 
by 



AP = P, . - P « = 1 AP 

actual ° 1 H i ideal 

where varies from 0.03 and 0.05. = 0.04 was selected as a 

representative average. 

Thus, the pressure drop becomes 



^actual ‘ ^ 






EQN (29) 



INLET AND EXIT GUIDE VANES 

The inlet guide vanes for the symmetrical blading were designed to 
provide the required flow direction, a^, to the first rotor stage 
at the design condition. The C-4 circular-arc profile was used. 

A solidity of 1.0 at the hub varying to about 0.85 at the tip was 
selected, with a chord varying linearly from hub to tip. The turn- 
ing angle variation with radius depends not only on the desired 
flow direction but also on the deviation angle of the flow leaving 
the guide vane. The method selected for determining this deviation 
angle is based on the two-dimensional cascade data of Dunavant (8) 
on a similar six percent maximum profile thickness guide vane, to- 
gether with the effect of axial velocity variation suggested by 
Soundranayagain (9) • A summary of the design data is given in 
Table 3. 
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The exit guide vanes for the symmetrical blading were designed to 
turn the flow leaving the last stator row back to the axial direc- 
tion, The radial variation of the inlet flow angles to the exit 
guide vanes is the same as the outlet flow angles produced by the 
inlet guide vanes. The maximum turning angle occurs at the outer- 
most or tip radius, . A solidity of approximately 1.2 was 
selected at the tip, varying to about 1.0 at the hub to provide a 
margin large enough to avoid flow separation. A linear variation 
of the blade chord from hub to tip was selected for the exit guide 
vane using the C-4 circular-arc profile sections. The aerodynamic 
design followed the method presented in Reference (1). A summary 
of the design data is given in Table 4. Note that the directional 
flow changes that occur in the exit guide vanes are equal and 
opposite to those that occur in the inlet guide vanes. Therefore, 
neglecting the minor total pressure drops occurring in the guide 
vanes, the magnitude of the pressure drop in the inlet guide vanes 
is equal to the magnitude of the pressure recovery in the exit 
guide vanes. 



DIFFUSOR 

The ideal pressure recovery in the diffusor at the exit of the com- 
pressor is given in terms of the inlet velocity head, 

4 - 4 V 8' - 1/2 V l 

and the diffusor area ratio, 

A • 

8 



23 



AP. 



** o 

1 ^ 1 \ 1 8 v 

D “ (P 9 “ V ideal " q 8 (1 ~ aT * 

ideal 9 



i z 

Note: Ag = A rgf and q g = q rgf (j) , therefore, 



AP 



P ideal 



[ ■ 0 ' 



1 ref 



f 



EQN (30) 



The actual static pressure rise can be defined in one of two 
ways. First, as a fraction of the ideal pressure rise, that is, 

ap d = n d ap d 

actual ideal 

where is the diffusor efficiency. Or, second, as a fraction 
of the maximum possible pressure recovery. The maximum possible 
is the recovery of the inlet dynamic head, or equivalently, when 
A approaches infinity. 

y 



AP 



D P 9 P 8 n D AP D. , n D q 8 n D q ref 9 ECJN 

actual ideal 

max 



where is the diffusor pressure recovery factor. Note that 
the two definitions are related to each other by the area ratio 



term, that is, 



n D = ( 1 



- («) 



) n • 



The latter method was chosen in order to use the results of Concanover, 
Kline, and Johnson (10) directly. 

A straight-sided, annular, conical diffusor with an area ratio of 
1.8 and a length of 2.4 feet was selected. From Reference (10), the 



recovery factor for the expected range of flow rates is approximately 
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0.56. (This corresponds to a diffusor efficiency, r = 0.81.) 
Two other diffusors with shaped contours and with approximately 
the same area ratio were examined and found to have recovery 
factors within ten percent of the selected diffusor. Because 
of the relative simplicity of manufacture, the straight-sided, 
annular, conical diffusor was selected. Note that the discharge 
pressure, p^ , is equal to the atmospheric pressure, P^ . 



OVERALL PRESSURE COMPATIBILITY 



A determination of the operating points for zero or no throttling 
can now be made with the components selected. The required com- 
pressor pressure rise, AP^ , is obtained by summation of the actual 
pressure changes (rises or drops) in the other components of the 
flow system. Note first that for essentially incompressible flow, 
the axial velocity components through the compressor remains con- 
stant. Thus, the increase in total pressure across the compressor 
is equal to the increase in static pressure, that is, 



AP = P - P = P ' - P ' = AP 
t t-i t.f 76c 

c 7 6 



Therefore, (see Figure 14) 



AP c = AP n - (AP e - AP T1 ) + AP th + AP t2 + AP X 



+ ap igv ap egv ap d 



EQN (32) 
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or 



P 7 ~ P 6 ^a ~ P 1 ^ 2 Pi ) + (P 2 " P 3 ) + Cp 3 - p 4 ) 

+ (p 4 ' - P 5 ') + Cp 5 ' - P 6 ’) - (p 8 ' - p/) - (p a " P 8 ') 



For zero throttling, Ap TR = 0. Also, Ap IGV = Ap^. 

Substitution of the previously obtained expressions for the compo- 
nent pressure changes and introducing the previously defined dimen- 
sionless pressure rise, it , into Equation (32) yields 

3 



2 Z tt q - = q c 
t^ ref n ref 



f* [ 



<1 - 5 „) - 2(1 - 1)(1 ' “ 2)2 

e e 



2 2 1 

+ 5 t2 — £ (1 - “ ) + 



(1 - tj (1 - (1 - “ 2) " %] 



EQN (33) 



or 



2 Z 



\ = ( F <s i' 3 ' 27 “ • “> e> V} f 



EQN (34) 



where F is the functional relationship given in explicit form in 
Equation (33) . 

Substitution of the previously selected values (summarized here) into 
Equation (33) 



’T 2 



0.04; 

0.04; 

0 . 02 ; 



Cj. = 0.04; 

n D = 0.56; 

a = 0 . 6 ; 



2 R 



= 3.0 



yields for 6 = 0.7 



F = 1.038 



EQN (35) 



for 8 = 1.0 



F = 0.506 



EQN (36) 
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These two expressions are plotted in Figure 15. Plotting the com- 
pressor performance on the same figure yields a graphical solution 
to Equation (34) . Also noted on Figure 15 are the design point and 
stall limit flow functions, 

Operation with one stage of symmetrical blading will be limited to 
(p ~ 0.95 which is higher than the design point jp = 0.825, but 

about five percent less than the predicted upper stall limit of 
<P ~ 1.03. Operation with two stages is possible throughout the 

complete flow range between the predicted stall limits with an inlet 
nozzle diameter ratio, 3 = 1.0; and for operation with three stages 
with a 3=0.7. 

Operation with one, two, or three stages of free-vortex or solid- 
body type blading is possible throughout the expected flow ranges, 
using an inlet nozzle with a diameter ratio, 3 = 0.7. 

Thus, the two inlet nozzle geometries selected will be adequate for 
the expected operating ranges. 



THROTTLING DEVICE 



From Figure 15, the difference between the compressor pressure 

2 

curves, 2 Z tt , and the sum of the losses curves, F <p , represents 
fc 3 



the pressure loss coefficient required by the throttling device. 
That is, 





n 

l 

i=l 




2 

2 Z tt - F <p 
t 3 J 



EQN (37) 
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n 



is the sum of n throttling elements in series, 



where E, . 

1-1 X TH 

each having a pressure loss coefficient of £ . . This required 

1 TH 

sum is plotted in Figure 16 as a function of the mass flow function, 
( p . The flow obstruction throttling elements selected are of 
two types: screens and per formated plates. The elements are circu- 

lar discs with the same inner diameter as the compressor inlet duct, 
and are welded to annular rings for rigidity. 

The pressure loss coefficients for both type elements were obtained 
from Reference (11). For screens, a correlation is found relating 
the pressure loss coefficient to the screen solidity, s , and the 
flow Reynolds Number based on the approach velocity and the screen 
wire diameter. The loss coefficient is given by 

C th - k«e> [ 1.3 S + (1 _ S g)2 ] 



where k(Re) is a function of Reynolds Number and varies between 0.9 
and 1.2 for the geometries selected. A summary of the screen con- 
figurations selected is given in the table below: 



Mesh Number, m (in. 


4 


10 


12 


12 


12 


Wire Diameter, d (in.) 

w 


.028 


.016 


.018 


.025 


.028 


Solidity, s 


.211 


.294 


.392 


.516 


.564 


Open Area, (1 - s) 


.789 


.706 


.608 


.484 


.436 


Loss Coefficient, L 

in 


.346 


.556 


.925 


1.807 


2.407 
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The material selected for the screen cloth was stainless steel. 

For perforated plates, it was determined that the feasible combina- 
tions of plate thickness, 1 , hole diameters, d , and solidities, 

n 

s , result in a range of V, between 0.2 and 1.0, and a range of 

H 

Reynolds Numbers (based on velocity through the holes, and the hole 

4 5 

diameter,) between 10 and 10 . Within this range of parameters, it 
has been found that the pressure loss coefficient for a desired plate 

solidity, and Re can be related to the pressure loss coefficient 

H 4 1 

for that solidity at reference values of Re = 5 x 10 , and = 0.4. 

H 

That is 

5 th f > ’ S TIi (=. *V </> 0 V K«e) ' K(j4) 

H H H 

where Re = 5 x 10^ and (^/ ) = 0.4. 

° d H ° 

A summary of the performated plate configurations selected is given in 
the table below: 



Hole Diameter, d (in.) 

H 


0.25 


0.625 


0.750 


Hole Spacing, (in.) 


0.375 


0.875 


1.0 


Plate Thickness, 1 (in.) 


0.1875 


0.1875 


0.1875 


Solidity, s 


0.60 


0.54 


0.49 


Loss Coefficient, 5 

(at Re ) TH 

o 


5.26 


4.37 


3.10 



The internal stresses and deformations of the screens and plates selec- 
ted were also investigated and found to be well within allowable limits. 
Similarly, the annular support rings were found to be well within allow- 
able limits to support the axial force exerted on the throttling elements. 
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FINAL ASSEMBLY 



Figures 17 through 21 are photographs of the assembled, redesigned 
compressor facility. Figure 17 is an overall view of that portion 
of the facility which is housed within the laboratory building. 
Figure 18 shows the inlet nozzle which is located outdoors. Figure 
19 shows the protective netting surrounding the inlet nozzle to 
minimize the effects of external wind gusts. Figure 20 shows the 
interchangeable throttling elements, one screen is shown and one 
perforated plate is shown partially withdrawn through the access 
door. Figure 21 shows the drive motor, pulley drive, and the 
torque meter installation. 
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Data 


at Mean 


Vavra Data 


Modified 


Velocity 






Radius, R 
m 


Ref (1), T .VII 


Triangle 


Results 


ROTOR 


INLET 


v :i- 


^al/a)R 

m 


0.8138 


0.8250 












0.3171 


0.3084 












0.6829 


0.6916 








< 




0.8734 


0.8807 












1.0624 


1.0765 








a i 




21.288 


20.496 








6 i 




40.0 


39.974 








AV* = 
u 


V a2/o)W* 

u 


0.38323 


0.38323 




ROTOR 


EXIT 


* 

V a2- 


^a2/o)R 

m 


0.8380 


0.8250 








* 














V u2 




0.7003 


0.6916 








<2 




0.2997 


0.3084 












1.0921 


1.0765 








* 

W 2 




0.8900 


0.8807 








a « 
2 




39.886 


39.974 








g 2 




19.68 


20.496 





TABLE 1. Comparison of Results of Modified Velocity 
Triangle with Those of Vavra (1) . 
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, original 
^1 

modified 


30.0 


32.0 


36.0 


40.0 


45.0 


50.0 


29.974 


31.974 


35.974 


39.974 


44.974 


49.974 


i' - i 


-10.0 


-8.0 


-4.0 


0.0 


+5.0 


+10.0 


. f 
1 


-11.828 


-9.828 


-5.828 


-1.828 


+3.172 


+8.172 


(d6/di) ref 


0.060 


0.062 


0.070 


0.078 


0.088 


0.097 


6 


9.79 


9.89 


11.11 


10.39 


10.83 


11.36 


» 

AB original 


10.92 


12.82 


16.60 


20.32 


24.88 


29.35 


ab'/ab ” 


0.537 


0.631 


0.817 


1.00 


1.224 


1.444 


Ag f modified 


10.46 


12.29 


15.91 


19.478 


23.84 


28.13 


Ci -i)/Ag mod. 


-0.492 


-0.394 


-0.197 


0.0 


0.246 


0.492 


t 

modified 


19.51 


19.68 


20.06 


20.496 


21.133 


21.85 


v* T = V* ? = (D 
al a2 7 


1.0520 


1.0020 


0.9094 


0.8250 


0.7284 


0.6392 


AV* 

u 


0.2339 


0.2670 


0.3280 


0.3832 


0.4462 


0.5048 


(w*) 2 


1.4749 


1.3951 


1.2627 


1.1589 


1.0601 


0.9878 


(v*) 2 


1.5001 


1.4155 


1.2731 


1.1589 


1.0467 


0.9080 


y'/Y 


2.21 


1.73 


1.14 


1.00 


1.17 


2.26 


Y 


0.0625 


0.049 


0.032 


0.028 


0.033 


0.064 


— t 

n t 


0.6025 


0.7421 


0.8753 


0.9144 


0.9221 


0.8798 


\3 


0.1345 


0.1891 


0.2740 


0.3345 


0.3927 


0.3927 


l'/L 


0.778 


0.846 


0.943 


1.00 


1.028 


0.946 



TABLE 2. Summary of Off-Design Calculations 
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TABLE 4. Summary of Exit Guide Vane Design Data. 
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Figure 1. Original Axial Flow Compressor Configuration 
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Figure 2. Redesigned Axial Flow Compressor Configuration 




Figure 3(a). Original Velocity Triangle (From Vavra (2)) 




Figure 3(b). Modified Velocity Triangle 
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Figure 4. Comparison of Modified Velocity Triangles 
at Design and Off- Design Conditions 
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Figure 5. Inter-relationship Between Geometric and 
Flow Angles for Compressor Rotor Blade 
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Figure 6- Normalized Modified Velocity Triangle 







Figure 7. Turning Angle as a function of Incidence Angle 
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LOSS COEFFICIENT RATIO 







Figure 8. Assumed Loss Coefficient Variation with Incidence Angle 




Figure 9. Position Locations within Compressor 
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Off-Design Pressure Ratio, Stage Efficiency, and 
Required Power for the Symmetrical Blading 
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Figure 10. 
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Figure 11. Off-Design Pressure Ratio, Stage Efficiency, and Power 

Required for the Free-Vortex and Solid-Body type Bladings 
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Figure 12. Inlet Nozzle Pressure Drop as a Function 
of Mass Flow Rate and Nozzle Throat Area 
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Figure 13. Inlet Nozzle Contour 
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Figure 14. Pressure Variation Through the Compressor Facility 
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Figure 15. Graphical Solution for Zero Throttling 
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Figure 16. Required Sum of Throttling Element Loss Coefficients 
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APPENDIX A 



PROCEDURE FOR OBTAINING THE MODIFIED VELOCITY TRIANGLE AT DESIGN POINT 
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